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ABSTRACT 


Hydraulic ram is the phenomenon associated with the dam- 
age caused by high speed projectile impacts with fuel tanks. 
It is divided into five phases, entry, shock, drag, cavity 
and exit. The entry, shock and drag phases were of concern 
їп this study. 

An analytical solution for the pressures developed in 
the shock phase was formulated based on bullet impact para- 
meters. 

Pressures developed in the drag phase were experimentally 
correlated with theoretical predictions. These pressures 


were proven to be dependent on front wall wave reflections. 
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I. INTRODUCTION 


The study of aircraft survivability has become increas- 
ingly important in recent years. This increase in impor- 
tance is due to several factors. Among these are the esca- 
lating cost and concomitant complexity of modern aircraft 
making the loss of an aircraft less acceptable. The in- 
creased complexity of modern aircraft also tends to make 
them more vulnerable to anti-aircraft weapons. At the same 
time anti-aircraft defensive capability is increasing; wit- 
ness the widespread introäuction of hand held infrared 
guided missiles of the Redeye class. These factors make 
survivability analysis and design of economic and military 
importance. One of the most promising areas for survivabili- 
ty analysis is the fuel system; in particular fuel cells. 
This is because the fuel cells occupy a large percentage of 
the total volume of the aircraft and are therefore most 
susceptible to damage. Damage to the fuel cell can result 
in fuel starvation, structural damage, fires or explosion; 
any of which can result in the loss of an aircraft. 

The interaction of the projectile and the fluid in the 
fuel cell produces damage that is many times greater than 
that caused by impact with an empty cell. One of the 
mechanisms that can cause damage is fire and/or explosion 
of the fluid. Hydraulic ram refers to the complicated 


dynamic interaction of the projectile, cell walls and the 
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fuel. These interactions produce internal pressure pulses 
which load the fuel cell structural members and produce 
damage. 

Hydraulic ram is usually divided into several phases 
that are associated with different major pressure pulse 
characteristics that occur uring the passage of the projec- 
tile through the fluid. The five phases of hydraulic ram 
are the entry, shock, drag, exit and cavity phases. The 
entry and exit phases occur when the projectile penetrates 
the respective fuel cell walls. The shock phase is associ- 
ated with the initial transfer of the projectile kinetic 
energy to the fluid in the cell. This energy transfer 
creates a strong shock wave. The drag phase occurs as the 
projectile continues through the fluid and transfers more 
energy to the fluid through viscous dissipation and pressure 
drag. The cavity phase occurs as the cavity, which was 
created in the drag phase, overexpands, is compressed, re- 
bounds and continues to oscillate until the energy is 
dissipated. The entry, shock and drag phases will always 
occur in order. The exit phase may not occur, and it may 
precede or overlap the cavity phase depending on the size of 
the tank, projectile impact velocity and other factors. 
Hydraulic ram research has been conducted by various investi- 
gators. This study was designed to extend Yurkovich's 


[Ref. 1] analytical solution for the shock phase fuel cell 
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pressure loadings and to verify experimentally the analyti- 


cal solution of drag phase pressures of Lundstrom [Ref. 2]. 
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II, DESCRIPTION OF HYDRAULIC RAM PHENOMENON 


In analyzing hydraulic ram all five phases must be con- 
sidered. The damage associated with hydraulic ram is depen- 
dent upon the anti-aircraft threat and the physical charac- 
teristics of the fuel cell. The threat determines the path, 
energy level and size of impacting projectiles. The fuel 
cell size, geometry, location and construction, in conjunc- 
tion with the threat, determines the dominant phase of the 
hydraulic ram phenomenon. For example, there may not be an 
exit phase for a low energy threat impacting a tank with 
sufficiently large dimensions along the projectile path. 

In any phase orf hydraulic ram the damage is related to 
the energy dissipated in that phase. This makes the deter- 
mination of the energy loss during the entry phase (wall 
penetration) essential. The energy dissipated during the 
entry phase is a function of wall parameters (thickness, 
yield strength) and projectile initial energy, obliquity of 
impact, material properties, and tumbling. Studies by Soper 
[Ref. 3] have shown that the percentage of energy loss for 


untumbled, normal entries can be approximated by: 


EZE 
ї ШВ _ 200 gy АТ, 1 (ТІ-1) 


° eee = Е; m М. 2 


1 
The studies also have shown that the penetration energy loss 


15 independent of projectile shape. If the projectile has 
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not tumbled and strikes normal to the wall, the damage area 
(A) approximately equals the cross sectional area of the 
projectile. 

If the projectile energy is sufficiently high, it will 
penetrate into the tank fluid and produce a strong hemis- 
pherical shock wave centered about the penetration point. 
This is the start of the shock phase. The shock phase is 
characterized by high pressures of short duration (microsecs) 
that are independent of fuel cell geometry. After the bullet 
enters the fluid, it has little effect on the shock wave 
strength or shape. Investigations have shown that shock 
phase pressures generated are a function both of the energy 
transferred to the fluid at impact and of the fluid proper- 
ties. The shock wave pressures are attenuated rapidly as 
the wave is propagated through the fluid. Therefore, the 
damage associated with the shock phase is generally limited 
to the entry wall. This damage 1s caused by the impulsive- 
ly generated pressure field combined with whatever damage 
(cracks, etc.) the initial penetration of the wall may have 
caused. 

An analytical formulation of the shock phase problem has 
been postulated by Yurkovich (Ref. 1]. In this formulation 
the shock formation process of the shock phase has been 
modeled as a point energy source release. The formulation 
assumes that the fuel cell walls are rigid, and that the 


shock radius is proportional to a constant power of time. 
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The assumption that the radius is proportional to a constant 
power of time has been supported by several investigators 
(Кес. 4, 5, 6]. 

The Rankine-Hugoniot equations are used to formulate the 
flow properties at the shock front. A fourth equation is 
provided by the Tait equation of state for isothermal com- 
pressible liquid. The fluid property changes at the shock 
front are assumed to be adiabatic. If the shock radius is 
known, a shock Mach number can be computed from energy con- 
Siderations. From the shock Mach number the shock front 
pressure and density can then be computed, using the above 
equations. Reference 2 postulates the initial radius is a 


function of energy deposited given by: 


D 
> ) 
Po C 27 


a = fF 1/3 (II-2) 


0 


where E, is the energy released in the shock phase, C is the 
speed of sound and 0g is the density in the undisturbed 
Fluid. 

The flow parameters behind the shock front are com- 
puted using the unsteady motion, one-dimensional equations, 
for conservation of mass, energy and momentum, again coupled 
with the state equation. To use these equations a power law 
density distribution was assumed. The radial velocity 


distribution was found by integrating the continuity equa- 


tion. Integration of the momentum equation then yields the 
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pressure. Integration of the energy density behind the shock 
front, knowing the shock Mach number and radius, would then 
yield the energy released during the shock phase. Since 

the energy released is assumed to be known, solving the 
energy integral for the shock Mach number is possible. Once 
the shock Mach number is known as a function of radius, the 
shock radius as a function of time can be found by integra- 
tion. It is then possible to determine the pressure, shock 
Mach number and radius, as a function of time, anywhere in 
the field. In this formulation, then, the energy released 
in the shock phase must be known. 

Since, in most cases, not all of the projectile energy 
has been released during the shock phase, another phenomenon 
of hydraulic ram becomes important. This phase is referred 
to as the drag phase. It is distinguished from the shock 
phase by the gradual acceleration of the fluid vice the 
impulsive acceleration characteristic of the shock phase. 
The energy transfer, therefore, is due to viscous dissipa- 
tion and form drag. Drag is normally characterized by 
defining a drag coefficient. The drag coefficient of a 
projectile is not a constant but varies with velocity, 
shape and tumbling behavior. If the projectile velocity is 
below the drag divergence Mach number, however, changes in 
velocity have only a small effect on drag coefficient size. 
The effect of shape and tumbling is then the dominant 


variable for estimating the drag coefficient for a 
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projectile. Experiments [Ref. 21 have indicated that soon 
after entry (within one foot) the projectiles will begin to 
tumble. Furthermore, because of the high stagnation 
pressures on the nose of the bullet, soft nosed projectiles 
may begin to deform. This phenomenon effectively changes 
the shape of the projectile and therefore the drag coeffi- 
cient. 

The pressure disturbance in the drag phase is propagated 
through the fluid at the fluid sonic velocity. Since the 
energy transfer is more gradual, the peak pressures are 
lower during this phase. But the pulses are of much longer 
duration (hundreds of microsecs). During this phase a cavi- 
ty is formed behind the projectile. The cavity is cylindri- 
cal in shape and its size and expansion rate are dependent 
on the rate of energy dissipation by the projectile. 

Lundstrom [Ref. 2] has formulated an analytical solution 
for the drag phase pressures. In this model a velocity de- 
cay factor B is defined as a function of the drag coefficient, 


mass, area, and density of the fluid: 
= = 
В = p_C, A/2m (11-3) 
Then the rate of energy dissipation with distance is given 


as: 


a= Se a) 


e 





о 


where B is assumed to be independent of velocity. The flow 


field properties are then solved for by using a potential 
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function that satisfies the wave equation. A line of sources 
is distributed along the projectile path with the strength 

of the sources being determined by an energy balance. The 
effect of wall restraint on cavity growth is not considered, 
so that the analysis will not yield useful pressures during 
cavity collapse. 

The pressure wave reflections at the fuel cell boundaries 
are assumed to be similar to those of a free surface. A sys- 
tem of negative image sources is utilized which result in 
zero pressure perturbation at the walls. Pressure waves re- 
flected from the cavity surfaces are not considered. 

The next phase to be considered is the cavity oscilla- 
tion phase. This is a complex phase and the least under- 
stood. Part of the complexity is due to the many reflected 
pressure waves that have had time to be of significance. 

The cavity is ae the drag phase as the momentum 

of the fluid becomes large enough to cause the fluid to 
separate from the projectile surface. A void is creäted 
which is filled with fluid vapors and air that enters through 
the wall penetration point. The cavity then grows until it 
is arrested by fluid pressures developed by the walls, and 
collapses. During the collapse the momentum of the fluid 
causes the entrapped gases and air to be compressed and the 
cavity begins to expand again. It is this alternating ex- 
pansion and contraction that produces low amplitude pressure 


pulses of long duration. The phase continues until all the 
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energy is absorbed by the fluid. High speed motion pictures 
[Ref. 7] have shown that the pressure oscillations can lead 
to eventual catastrophic failure of walls already weakened 
by the drag and exit phases. 

The exit phase will occur if the energy transferred from 
the projectile during the entry, shock and drag phases has 
not dissipated all the kinetic energy of the projectile by 
the time the exit wall is reached. Pressure loadings due to 
projectile proximity to the wall, and fluid motion produced 
in the earlier phases, can cause catastrophic damage to the 
exit wall. If this occurs, the latter phase (cavity oscilla- 
tion) causes little further damage but does pump fuel from 
the tank, causing fire hazard. 

From the above discussion it is apparent that each of 
the different phases of hydraulic ram produces damage to the 
fuel tank in a eT manner. The entry phase results in 
a hole in the tank entry vall. The damage from the shock 
phase is usually limited to the entry wall and can lead to 
catastrophic failure. The drag phase results in lower 
pressure amplitudes, but of longer duration, which can cause 
high levels of stress on all fuel cell walls. The cavity 
phase has even lower pressure amplitude but longer pulse 
widths and is oscillating in nature. This can result in 
damage to the already weakened fuel cell walls. The exit 
phase causes damage to the exit wall, usually before the 


cavity phase. It should also be apparent that impact energy 
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and the physical characteristics of the fuel cell determine 


which phase is dominant. 
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ESI HYDRAULIC. RAM PRESSURE EXPERIMENT 


The Naval Postgraduate School Ballistic Range was used 
to conduct the experimental part of this study. The basic 
components of the range were a rifle, velocity measuring 
equipment, simulated fuel tank, and pressure measuring and 
recording equipment. Figure III-l is a schematic of the 
basic components and Figure III-2 is a down range view. 

The rifle was a .222 caliber Remington. Cartridges were 
loaded to obtain the 7,493 in.-lb. energy level with a 45 
grain spitzer ogive nose bullet. The velocity of the bullet 
was measured by using the first three Avtron No. A914T33 
chronograph screens to trigger two Monsanto 101B counters. 
Two counters were used to increase the accuracy of velocity 
measurements and for redundancy. 

The simulated fuel tank was constructed of welded .160- 
inch 7075-T6 aluminum plate on the sides, back and bottom 
walls. The front wall was .05-inch 7075-T6 aluminum. А опе- 
inch hole was drilled in the front wall for bullet entry. 
The hole allowed calculation of precise impact energy. The 
top of the tank was left open. The side and bottom walls of 
the tank were welded to a frame of .5-inch thick, three-inch 
aluminuin angle sections for added strength. The entry and 
exit wall were bolted to the frame. The inside dimension of 
the tank was 18x18x18-inches. Figure III-3 shows the tank. 

A Kistler 603A quartz pressure transducer was used to 


measure the hydraulic ram pressure pulse at several points 
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internal to the tank. The transducer was mounted in a 

hollow 1.0-inch stainless steel tube. The transducer was 
held in the tube by a nylon ring. The tube was then sealed 
with RTV. In this manner the transducer and cable were kept 
dry. The tube had a 45° angle at the tip and was mounted on 
the tank at 50° from the vertical to obtain pressure measure- 
ments along a radial line. Figure III-4 shows the transducer 
and probe assembly. The probe was then secured to the top 

of the tank so that it could be positioned for measurements 
at 2,4,6,8 and 10 inches along a radial line from the impact 
Point. 

The charge signal from the pressure transducer was fed 
into a Kistler Model 504E charge amplifier with a 545A14 
filter and displayed on an oscilloscope. A Polaroid camera 
mounted on the oscilloscope recorded the pressure trace. 

The oscilloscope was triggered by the fourth chronograph 
Screen mounted one inch in front of the tank wall. 

Several firings were made at each distance for consis- 

tency and accuracy. Typical results for each distance are 


shown in Figures III-5, 6, 7. 
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Radius = 2 in 
Scale: 1000 psi/cm, 50 microsecs/cm 





Radius = 4 in 
Scale: 500 psi/cm, 50 microsecs/cm 


Figure III-5. Pressure vs. Time, Eg = 7,493 in-lb Light 
Wall. 
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Figure III-6. Pressure vs. Time, =" 7,493 in-1b Light 
Wall. Scale: 500 P ui 50 microsecs/cm. 
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Radius = 10 in 


Scale: 500 psi/cm, 50 microsecs/cm. 


Figure III-7. Pressure vs. Time, E, = 7,493 in-lb Light 
Wall. н 
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IV. ANALYSIS 


The first segment of this study was to improve the com- 
puter code for shock phase pressures deveioped by Kappel 
(Кес. 4]. This code was based on the formulation of 
Yurkovich discussed in Section II. 

To account for the energy dissipated by the projectile 
in penetrating the entry wall the formulation of Soper was 
added to the computer code. 

A method for computing the energy deposited during the 
shock phase was needed as previous formulations assumed 
this value was known. Studies [Ref. 4, 5] have indicated 
that the energy released in the shock phase is a constant 
percentage of the impact energy. This percentage is a func- 
tion of the impact Mach number with respect to the fluid 
speed of sound. Plots of experimentally observed changes 
in projectile velocity during the characteristic shock 
formation with projectile impact were made by Power [Ref. 8]. 
This plot is reproduced as Figure IV-i. In addition, the 
corresponding shock phase energy dissipation percentage of 
the impact energy is shown. It can be seen that impact Mach 
number with respect to the tank fluid is critical in deter- 
Mining the shock phase energy input. For Mach numbers less 
than sonic the input percentage is approximately 10 percent 
and for moderate supersonic speeds the energy input percent- 


age is greater than 80 percent. This provided an empirical 
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procedure for determining shock phase energy input as a func- 
tion of bullet impact parameters. 

The shock front pressures are quickly attenuated and be- 
come acoustic. Figure IV-2 shows the radius and time for 
the shock wave to become acoustic versus input energy for 
both water and fuel. To calculate pressures after the 
acoustic speed was reached, a simple expression based on 
acoustic wave theory was used. The shape of the wave was 
assumed to remain constant and the pressure amplitude varied 
as the inverse of the radius. 

Since the Yurkovich theory postulates the energy release 
aS a point source, the model can be extended to an explosion 
in the fuel cell. This type of event is characteristic of 
an H.E. round impact. This was accomplished by integrating 
the energy equation over a spherical volume. In this case, 
however, the shock phase energy that is released by the explo- 
Sion must be known. 

Figures IV-3, 4 and 5 show the predicted pressure dis- 
tribution for .90 caliber, 12.7mm and 14.5mm AP projectiles 
impacting JP-5 fuel. The plots illustrate the very high 
pressures РАЙ and their rapid decay. Figures IV-6 
through ll show the predicted shock phase pressure versus 
time at a constant radius for a .222 caliber projectile im- 
pacting water. Also shown in these figures are measured 
pressures from reference 9 and analytical drag phase 


pressures to be discussed. From these figures it can be 
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seen that the shock phase pressures could not be measured. 
This is due to the very rapid rise time of the shock phase 
pressures. The pressure transducers and associated measuring 
equipment were unable to follow the signal. This phenomenon 
was also reported by Williams [Ref. 10]. 

Experiments were undertaken to verify the drag phase 
pressures predicted by Lundstrom. Holm in reference 9 re- 
ported that the drag phase pressures couid not be verified. 
In that experiment the front wall of the fuel tank was a 
heavy, .5-inch, mild steel plate.  Lundstrom's formulation 
of the problem treats the walls as free surfaces. This is 
done by using a set of negative image sources so that the 
pressure at the wall is zero. 

The degree to which a wall may be considered rigid or 
free is characterized by the expression [Ref. 2]: 

а= рс С (IV-1) 
where 0 is the time constant. A rigid surface corresponds 
to a=0, and a free surface to О-о. For a .5-inch thick 
steel plate a=1.5, while for a .05-inch thick aluminum plate 

=44.5. The heavy steel plate is then an approximation to 
a rigid surface while the .05-inch aluminum plate is an 
approximation to a free surface. 

A set of image sources that can approximate the rigid 
surface would be a positive set. Using this set of positive 


images on the front wall for the pressures reported in 
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reference 9 resulted in very close agreement. This is 
illustrated in Figures IV-6 through 11 for two energy levels. 
The peak pressure versus radius is depicted in Figure IV-12. 
Close agreement is shown at medium and larger radii. The 
theory does underestimate the pressures at small radii how- 
ever. 

The experiment was repeated using a .05-inch aluminum 
front wall to approximate a free surface. Plots of the 
pressure distribution versus time at a constant radius are 
shown in Figures IV-13 and 14. These figures also show 
close agreement between the experiment and theory when a 
free surface is assumed. 

A plot of peak pressure versus radius compared with 
theoretical values is shown in Figure IV-15. Again, the 
theory underestimates the peak pressure for small radii. 

All of the theoretical pressures were computed using a 
drag coefficient of one. This coefficient was determined 
from the experimental data of reference 3. The theoretical 
formulation allows a variation of drag coefficient with 
time. This would account for any tumbling of the projectile. 
No measurements of the tumbling behavior were possible, so 
a constant Co was assumed. Photographs of the bullet 
trajectory would be required to determine tumbling distances. 

The pulse shape discrepancy was due to the assumption 
that C was a constant. If detailed measurements are made, 


the pulse shape will be more accurately predicted. 
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Reference ll reports that a data base was established for the 
tumbling behavior of 14.5mm and 12.7mm projectiles. From 


these data good correlation of pressures was possible. 
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47 





Pressure (psi ) 


Pressure (psi) 


200 


IOO 


200 


IOO 





[30 


ISO 


XXX Experimental 
—_— Calculated Drag Phase 


Radius = 8 in 


ISO (70 I90 210 230 
Time (p sec) 


Radius = IO in 


170 [90 210 230 
Time (y sec) 


FIGURE IV-14 PRESSURE COMPARISON 


Eo =7,493 IN.LB., LIGHT WALL 


48 





TIVM LH9N YOJ SNIQaVY4 “SA 3Ч4055344 ЯУ34 “St-AI 341913 


(seuoui) зтром 
OO! 06 O8 07 09 OS 0% 0% оа 


OO OI 


0002 


O O 


quuiC6p'4 soy бола рофојпојо) — 
аш бв‘ IpjueuiJjedx3 О 


(1$9) олп55оја 


49 





V. CONCLUSIONS 


The first purpose of this study was to extend the com- 
puter code solution for shock phase pressures. This was ac- 
complished by adding solutions for the entry phase and for 
the shock phase energy release. The attached code computes 
the shock phase pressures based on a knowledge of bullet im- 
pact parameters only. It will also compute the shock phase 
pressures for a point explosion in the fluid. 

Attempts to verify the shock phase pressures ee 
En were unsuccessful. This was due to the inadequate re- 
sponse time of the pressure measuring equipment. Until 
pressure transducers with extremely rapid rise times are 
available, shock phase pressures cannot be experimentally 
determined. 

The second part of this study experimentally confirmed 
the analytical solution for the pressures produced in the 
drag phase, postulated by Lundstrom. In accomplishing this, 
the importance of entry wall reflections on tank internal 
pressures was confirmed. By using both a heavy and light 
wall in the experiment, and matching these with appropriate 
image sources in the analytical solution, the pressures were 
correlated. 

The use of a constant drag coefficient in the analytical 
solution does distort the pressure wave shape. This is 
because the use of a constant drag coefficient is not con- 
Sistent with the tumbling behavior of the bullet. At 


50 





smaller radii the experimental pressures are higher than 
predicted. This is due to the interaction of the shock 
phase pressures with the drag phase. That is, the pressure 
transducers are partially responding to the impulsive 
pressure of the shock phase. 

The results of this study show that the first three 
phases of hydraulic ram are clearly understood and can be 
predicted analytically. This should be of some assistance 
to those responsible for increasing the survivability of 


aircraft subject to light or moderate damage. 
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ООС ОС 


MOOD 


= 


ЕТО) ВОО, 109,90 


SUBROUTINE SUBENLIFLAS,CO,RHOO,XN,VSE,ET) 


THIS SUBROUTINE COMPUTES THE ENERGY BEHIND THE SHOCK 
FRONT FOR A GIVENSMACH NO. THE EXTERNAL INTEGRATION 
SUBROUTINE OG1I IS A GAUSSIAN 10 POINT QUADRATURE. 
EXTERNAL FCT 

COMMON BQ 

C=2.0 

F CIFLAGIEQ.0) C=1.0 _ 

САШЕ SHOCK [АМ УЕ, РЪТЕ, Ре) 

A=3.1416* (RHO0/32.2)*[C0*C0/144+.0)*(VSE**2.0)*PSTE 
*FE*FE*2,0 

BQ=3.*PSIE-1 ° 

CALL ОБО 00.0.1. 0.ЕС Т, ЕЕ) 

ЕТ= А* ЕЕС 

RETURN 

END 

FUNCTION FCT(X) 

THIS FUNCTION SJ3PROGRAM GENERATES FUNCTIONS FOR THE 
INTEGRATION ROUTINES. 

COMMON BQ 

met—( te tx**2 . ) +X *BQ 

RETURN 

END 
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